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FOREWORD 

This  interim  report  was  prepared  by  Dr.  Gdalia  Kleinstein,  Research 
Associate  Professor  of  Aeronautics  and  Astronautics  at  New  York  University, 
and  presents  research  carried  out  under  a  National  Science  Foundation 
Institutional  Grant,  "Statistical  Methods  and  Their  Application  to  Fluid 
Dynamics  Problems",  and  Contract  No.  AF33  (615)-2215,  "Boundary  Layer 
Characteristics  for  Hypersonic  Flow  in  the  Presence  of  Mass  Addition", 
Project  No.  7064. 
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ABSTRACT 


Based  on  the  logarithmic  velocity  distribution  at  distances  not  so  far 
from  the  wall  and  a  turbulent  shear  variation  proportional  to  the  cubic  power 
of  the  wall  distance  in  the  immediate  proximity  of  the  wall,  an  analytical  model 
for  an  eddy  viscosity  throughout  the  wall  region  is  derived  leading  to  a  closed 
form  expression  for  the  velocity  distribution.  An  extension  of  the  above 
analysis  to  include  shear  variation  results  ir.  a  generalized  law  of  the  wall 
from  which  the  effect  of  Reynolds  number  on  the  velocity  distribution  is 
derived.  While  many  applications  of  the  generalized  solution  are  envisaged, 
as  an  example,  the  velocity  distribution  for  the  case  of  injection  or  suction 
is  considered.  The  theoretical  result  is  shown  to  be  in  agreement  with 
available  experimental  data. 
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NOMENCLATURE 


mean  flow  velocity  (streamwise  component) 
friction  velocity  (local)  =  |\(y)/pj^ 
wall  friction  velocity  u*  =  (u  )^  =  j\(0)/pj^ 
dimensionless  mean  flow  velocity  ur  «  u/u 

4-  ^ 

dimensionless  friction  velocity  u  =  u  /u 

T  T 

eddy  viscosity 

dimensionless  eddy  viscosity  t/v 
normal  coordinate 

4 

dimensionless  normal  coordinate  y  =  yu  /v 

4. 

dimensionless  injection  velocity  v  =  v  /u 


INTRODUCTION 

The  validity  of  the  logarithmic  nature  of  the  law  of  the  wall,  which 
may  be  derived  on  the  basis  of  substantially  different  underlying  assumptions 
such  as  dimensional  analysis  as  well  as  the  Malkus  theory*-,  has  been  satis¬ 
factorily  verified  by  a  large  volume  of  experimental  data.  Since  the  law  of 
the  wall  as  such  does  not  hold  through  the  transition  region  (buffer  zone)  and 
obviously  not  in  the  laminar  sublayer,  it  has  been  desired  to  obtain  a 
formulation  which  will  connect  continuously  and  smoothly  the  laminar  sublayer 

velocity  distribution  to  the  law  of  tne  all.  Most  successful  has  been 
2 

Reichardt  who  from  considerations  of  the  behavior  of  the  velocity  fluctuation? 

in  the  immediate  vicinity  of  the  wall  deduced  that  the  eddy  viscosity  must  grow 

at  least  as  the  cubic  power  of  y.  Accordingly,  by  choosing  a  function  which 

satisfies  the  y3  behavior  near  the  wail;  and,  asymptotically  approaches  the 

required  linear  dependence  on  y,  he  was  able  to  construct  an  eddy  viscosity 

model  which  in  turn  led  to  a  velocity  distribution  in  agreement  with 

experimental  data  from  y  =  0  to  the  law  of  the  wall  region.  Other  formulas 

34 

for  the  eddy  viscosity  have  been  suggested  by  Deissler  and  Van  Driest 
which,  .although  successfully  predict  the  velocity  distribution,  raise  the 
objection  that  they  do  not  satisfy  the  y3  requirement  of  Reichardt's  theoretical 
result. 

In  this  paper  it  is  shown  that  by  taking  the  velocity  as  the  independent 
variable,  an  eddy  viscosity  model  can  be  derived  from  the  law  of  the  wall  and 
Reichardt's  y3  decay  requirement  without  introducing  a  new  arbitrary  function. 
The  resulting  eddy  viscosity  model  yields  an  integrable  closed  form  expression 
for  the  velocity  distribution  as  compared  with  the  models  of  Reichardt,  Deissler 
and  Van  Driest  which  require  a  numerical  integration. 

By  employing  a  simple  transformation  the  above  formulation  is  further 
extended  to  include  the  effect  of  shear  variation.  Consequently,  a  generalized 
form  for  the  law  of  the  wall  is  obtained;  and  in  particular  it  is  shown  that 


1 


for  regions  where  the  eddy  viscosity  is  much  larger  than  the  kinematic 
viscosity  the  law  of  the  wall  has  a  universal  nature  independent  of  the 
shear  distribution.  The  application  of  this  result  to  pipe  flow  provides  an 
explanation  for  the  effect  of  Reynolds  number  on  the  velocity  distribution 
which  has  been  observed  by  Hinze^.  As  an  additional  example,  the  generalized 
formula  has  been  applied  to  the  problem  of  wall  injection  and  suction  uh°re 
good  agreement  with  experimental  data  has  been  demonstrated  in  regions  close 
to  the  wall  including  the  buffer  laye,.. 


ANALYSIS 


1.  Introduction 

The  total  shear  associated  with  the  transport  of  momentum  in  the  direction 
normal  to  the  wall  is  given  by 


du  — r-  / 

^‘PVU 


(1) 


If  a  concept  of  an  eddy  viscosity  is  adopted  the  turbulent  contribution  may 
be  expressed  by  the  relation 


-pv  u 


du 

P£¥ 


(2) 


where  Eq.  (2)  is  taken  as  a  definition  of  the  eddy  viscosity.  In  terms  of  the 
nondimens ional  quantities  y  =  (pu  y)/p,  u  “  u/u  ,  t  =  pt/p.;  where  u  *  (Tw/p) 
Eqs.  (1)  and  (2)  may  be  combined  in  the  form 


T  _  ,  +xdu 

jr  *=  (l+€  ) — ~ 

T  ay+ 


(3) 


w 


or 


(»!>=-  (i+.+)^r 

T  5y+ 


(3') 


4*  ★  \ 

where  u  =  u  /u  =  [t/t  is  the  nondimensionalized  friction  velocity. 

T  T  W 

Transforming  the  variables  from  (x,y)  to  (x,u)  yields  Eq.  (3  )  in  the 
desired  form  i.e. 


€+  +  1  =  (u+ 
T 


>=  8*! 


(4) 


du 


By  dimensional  analysis  considerations,  Prandtl  in  his  momentum  transport 
theory**  proposed  an  eddy  viscosity  model  which  in  the  present  notation  takes 
the  form 


*  k?(y  r  4 

dy 

where  partial  derivatives  have  been  replaced  by  an  ordinary  derivative  to 
indicate  a  dependence  on  a  single  variable.  Eliminating  du+/dy4  from  Eq. 

(5)  by  Eq.  (3  ;)  where  c+  is  taken  much  larger  than  unity  yields  an  equivalent 
model  in  the  form, 
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(6) 


\ 


4*  4*  4* 

«  *  h  v  (6) 

The  model  given  by  Kq.  (6)  has  been  previously  used  by  Ferrari^  in 
several  problems  and  specifically  in  problems  related  to  turbulent  boundary 
layer  separation.  Ecs.  (4)  and  (6)  in  conjunction  with  the  requirement  of 
(y+)3  decay  at  the  wall  constitute  the  basis  for  the  derivation  of  the  complete 
velocity  distribution  and  eddy  viscosity  formulation  given  below. 

2.  Constant  Shear  Solution 

If  the  shear  is  assumed  constant,  u*  *  1,  and  Eqs.  (4)  and  (6)  become 


and 


+  dy~* 

C  K  —  * 

du+ 


+  .  + 
<  m  hy 


(4') 


In  Eq.  (4  ')  t+  is  taken  much  larger  than  unity  and  a  dependence  on  a 
single  variable  is  assumed. 

Eliminating  between  Eqs.  (4 ')  and  (6 ,)  and  integrating  yields,  the 
logarithmic  law 

u+  =  -p  4ny+  +  ka  *  ^  inkay*  (7) 

N. 

where  k^  =  0.4  is  the  von  Karman  universal  constant,  and  ka  is  taken  equal 
to  7.7  which  is  obtained  from  ka*  =  5.1. 

Using  u+  as  the  independent  variable  Eq.  (7)  is  written  as 


+  1  ,,  +* 
y  =  j^expC^u  ) 

and  correspondingly  the  expression  for  £+  becomes 

€+  =  ^exp(l^u+) 

In  order  to  extend  the  above  solution  to  the  wall  (i.e.  y+  =  0)  the 


(8) 


(9) 


behavior  of  pv  u  and  u  must  be  specified  in  this  region.  For  the  immediate 
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vicinity  of  the  wall> Reichardt  has  shown  that  the  turbulent  shear  varies 
proportionally  to  at  least  the  cubic  power  of  the  distance  from  the  wall,  and 
since  for  the  same  region  the  mean  velocity  u  varies  linearly  with  the  same 


dimension  it  follows  from  the  definition  -pv  u  -  p«du/dy  that 


3  “3 

p€  ~  y  ~  u 


(10) 


t 


4 


i 


Now  including  a  correction  term  in  Eq,  (9),  t  may  be  written  in  the 


f"*  =  ^[exp^  u+)  -4  f  (kx u+)] 
where  f  must  satisfy  the  asymptotic  condition 


Urn  +  + 

y+  at  (kx  u  )exp(-k1u  )  =  0 

and,  in  accordance  with  Eq.  (10).  the  boundary  conditions 


(12') 


f(0)  =  f'(0)  =  f" (0)  «  -1  (12") 

While  Eq.  (12)  obviously  does  not  define  f  uniquely,  an  exponential 
representation  for  f  in  the  form  cx  exp(c3k1u+)  where  cx  and  Cg  are  arbitrary 
constants  is  clearly  excluded  by  Eqs.  (12')  and  (12").  Taking  then  the  highest 
degree  polynomial  which  satisfies  both  conditions  and  at  the  same  time  es 
not  introduce  new  arbitrary  constants  yields, 


f(*qu+)  ~[l  4-  (kxU+)  4  |(klU+)2] 


Accordingly,  the  modified  expression  for  t  becomes, 

«■*  =  ^-{Cexp (kxuH)3  -  Cl  4  (kxu+)  +  |(kxu+)2]}  (14) 

Inserting  Eq.  (14)  into  Eq.  (4)  with  u*  =  1,  and  integrating  subject  to  the 
boundary  condition  y+  =  0  yields  the  velocity  distribution 

y+  =  u+  +  -j^{[exp(k1u+)]  -  [1  +  (l^u*)  +  -^(1^  u-4")32  4  •|(k1u+)3]j-  (15) 

A  plot  of  Eq.  (15)  is  given  in  Fig.  (1)  where  the  classical  law  of  the 
wall  and  the  laminar  velocity  distribution  for  the  sublayer  are  included  for 
comparison.  Fig.  (2)  is  a  reproduction  of  Fig.  (3)  of  Lindgren's  recent 

g 

report  on  experiments  in  turbulent  pipe  flow  of  distilled  water  .  The  solid 
line  is  a  plot  of  Eq.  (15)  using  the  logj^y+  representation.  The  agreement 
with  these  measurements  is  surprisingly  good  within  the  range  0  <  u+  <  18. 

The  deviation  beyond  this  region,  although  small  in  the  beginning  becomes  more 
and  more  significant  as  the  wake  region  is  penetrated. 

A  comparison  between  as  given  by  Eq.  (14)  and  the  measurement  of 


as  tne  variable.  Taking  U  as  the  independent  variable, 


and 


y+  ' 


t+  =  u+  ^-exptkiU) 
T  Kg 


(21> 


(22) 


Now,  since  independent  of  the  nature  of  the  shear  distribution  there  always 
exists  a  narrow  region  close  enough  to  the  wall  where  U  may  be  taken  equal  to 
u+,  it  follows  from  Reichardt  condition  that 

«+~  (y+)3  ~  (U+)3  ~  (U)3 

and  consequently  by  the  identical  arguments  of  the  previous  section. 


Integrating  the  equation 


+  =  ytexpO^U)]  -  [l  +  (kjU)  +  |(kiU)a]} 


(23) 


+  ,  +  dy+ 

c  +  1  =  u  — rjr 
T  du 

subject  to  the  boundary  condition  y+  =  0  at  U  =  0  yields  the  generalized  law 

of  the  wall  in  the  form, 

U 

Y  *  J  +  HCe*P(M>]  -  [1  +  (krU)  +  IohU)2  +  Ic^U)6]}  (24) 

o  u 

T 

It  shoule  be  noted  here  that  there  is  a  fundamental  difference  between 
Eq.  (24)  and  Eq.  (15).  While  Eq.  (15)  represents  y  as  a  function  of  u+  only 
Eq.  (2*+)  nay  depend  on  the  Reynolds  number  as  for  example  in  a  pipe  where  the 
shear  is  given  by 


i_  =  i  _ 

T  d 

w 


(25) 


Now,  if  Eq.  (24)  is  to  predict  the  correct  velocity  distribution  for  a 

pipe  it  remains  to  explain  the  apparently  good  agreement  of  Eq.  (15)  with 

experimental  data  where  the  shear  was  assumed  constant.  Since  in  pipe  flow 

seventy  to  eighty  percent  of  the  centeriine  velocity  is  attained  within  a 

region  close  to  the  wall,  corresponding  to  2y/d  of  the  order  of  fifteen  to 

twenty  percent,  and,  since  beyond  this  region  the  flow  field  is  definitely 

within  the  region  of  the  law  of  the  wake  it  is  required  (1)  to  show 

that  E  q  .  (24)  and  Eq.  (15)  do  coincide  within  the  range  0  <  2y/a  <  0.20 
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to  a  close  approximation  independent  of  the  shear  variation  and  (2)  to  determine 
the  Reynolds  number  dependence  of  the  law  of  the  wall  in  a  pipe. 

In  terms  of  the  nondimencionalized  variables,  the  friction  velocity  for 


the  pipe  is  given  by 


u  -  a  -  '•  a  -  ^ 

T  c  d 


I  which  for  small  values  of  the  ratio  2y/d  may  also  be  written  as, 

u+  =  1  -  y  /d  4  0[  (y  /d  )2  ]  (5 

T 

+  +  + 

For  the  zero  order  solution  of  (24),  designated  by  yQ,  let  u^  =  (u  )q  ~ 
Correspondingly  U  *=  u+,  the  integral  term  reduces  to  U  and  therefore  to  u  ; 
and,  thus  the  velocity  distribution  of  Eq.  (15)  is  recovered.  For  the  first 


(26') 


order  solution  let 


4  /  +N  _  1  0 

UT  “  <u  >i  "  1  ’  T 
T  t  d4 


Hence, 


V* 

P  du'  P  /,  yo,  /  4 

U  -  J  =  j  (1  4  — ^-)du  =  u  4 


4  (d+)~1p|(u't)~  4  j^-|[exp(k1u+)j  "  L1  4  (ki u+)  4  y(kiu+>s  4  ^/kiu  )3  4 

4 1-/MV]}]  <27> 

and  taking  into  account  that  y+(u4)  =  y"|"(U), 

U  U  4  ° 

i  =  j1—  =  fa  4  )du'  =  u  4 

-  J  U  J  ,+ 

o  t  o  a 

4  (d-1)"1^2  4  |[exp(k1U)]  -  [1  4  (kxU)  4  |(kxU)3  4  |(kxU)3  4  ^(M)4]}} 

1  (28) 

With  U  as  a  known  function  of  u*  and  I  resolved  in  terms  of  U,  Eqs.  (27)  and 

(28)  in  conjunction  with  Eq.  (27)  completely  describe  the  dependence  of  the 

4  * 

velocity  distribution  on  the  Reynolds  number  d  =  u  d/v. 

As  an  immediate  consequence  of  this  solution  it  is  apparent  that  for 

■f  4  4 

any  fixed  value  of  u  ,  yx  approaches  the  constant  shear  solution  yQ  as  the 

Reynolds  number  d+  approaches  infinity.  This  result  may  also  be  directly 
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I 

f  5 

* 

inferred  from  Eq.  (20)  since  a  fixed  u+  implies  a  fixed  y*. 

In  order  to  assess  the  maximum  deviation  between  y  and  yQ  it  is  necessary 
to  consider  the  asymptotic  region  where  independent  of  how  large  d+  is,  the 
ratio  2y^/d  may  asrume  values  up  to  0.20.  For  the  asymptotic  region  Eq.  (27) 
and  (28)  reduce  respectively  to 

U  =  u+  -f  ~  ^  exp  (kx  u  )  (27 f ) 

and 

I  “  u  +  ~  A-  expO^U)  (28') 

(j"1"  K1  *3 

Substituting  Eq.  (28/)  in  Eq.  (24)  yields 

y4  -  u  +  (r$r  4  l)^-exp(k1U)  (29) 

d  kj^  ** 

where  for  the  time  being  quantities  of  the  order  of  (d+)  are  retained  with 
respect  to  U. 

Using  Eq.  (27)  to  eliminate  U  from  Eq.  (29)  yields 


4-  4  1  1  ..  4-.  ,  1  ,.l  4-.  rl  1 _  .  4.-, 

y  =  u  4  —  l^exp(k1u  )  4-  (-^—4  l)-exp(k1u  )exp[—  ^expO^ u  )] 

d  d  k^  d 

Now,  since  asymptotically  y+  =  7-exp  (kiU+)  the  above  equations  may  be 

O  42 


written  as 


4-  4-  1  'o 


+  7os 


y 1  ■”  +  *7+<a  +  °  ? 

'  “+  +  k  ^ +  +  1,y°(1  +  S H  0 

1  d  +d  ^  £  d 

-  <u+ + ,+) ♦ ♦  y ♦ c[<j)2] 

But  since  y+  >  >  u+  and  y  ^  c^e  asymptotic  value  of  y*  becomes 


*  =  >&' +  7f] +  { ("r>aj 


While  as  far  as  y i  is  concerned, this  deviation  may  be  as  large  as  ten  percent, 

as  far  as  the  velocity  distribution  is  concerned  where  uT  is  taken  proportional 

4-  4-4- 

to  tny  ,  the  deviation  Au  /u  given  by 


i 

t  * 


\ 


Au+  W^ny*  y^/d+ 

—  <  «=  —  **  ■  — 

u  4ny  <tny 

o  Jo 

is  approximately  within  one  percent. 

A  different  interpretation  of  this  result  may  be  given  as  follows.  Let 

+  +  + 
y  be  taken  as  a  function  of  d  then  at  a  fixed  value  of  y 

J  o 

ys  Cl  Cs  C]  cfe 

Thus,  for  >  d^  *»  y4  >  y t •  Furthermore .since  at  increasing  values  of  y4 
the  difference  between  y^  and  yt  increases.it  follows  that  as  the  Reynolds 
number  increases  the  slope  of  the  straight  line  in  the  u  ~  -tny  representation 
increases  and  correspondingly  the  intersection  with  the  y4  £  1  line  occurs  at 
a  lower  value  of  the  u+  coordinate.  In  terms  of  the  notation 

u"  =  ^ny+  +  ka 

it  means  that  1/kx  increases  while  ka7  decreases  with  an  increasing  Reynolds 
number.  This  dependence  on  the  Reynolds  number  has  been  shown  to  exist  in 
Nikuradse's  data  by  Hinze3. 

In  order  to  see  the  effect  of  the  Reynolds  number  on  the  laminar  sublayer 
a  series  expansion  of  the  solution  in  this  neighborhood  is  considered. 
Accordingly 

y+  =  u+  4 1  k(kiu4)4  4  <d+)_1t(u+)2  +  ito(l^u+)5  4---]  (31) 

By  comparing  this  result  to  the  strictly  laminar  distribution 
u  =  y  -  (d  )  y 

which  for  large  Reynolds  numbers  and  y/d  <  1/2  may  be  written  in  the  form, 


4  4  ,,4.-1  4a 

y  =  u  4  (d  )  u  4  0L(d  ) 


the  extent  of  the  laminar  sublayer  may  be  obtained  from  the  condition 


fc  24 U+>'  <  <  tfWf 
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Taking  as  the  definition  of  the  edge  of  the  laminar  sublayer  the  point  where 

-°  l(d4>'l<“s.L.)2  <32) 

yields,  after  introducing  the  numerical  constants 

Ug  L  ~  27 (d+)^ 

and  by  substitution  into  Eq.  (31) 

+  u*d 

yS.L.  -  27 


yS.L.  07  .u* dv -3/2  (33) 

Hence.it  follows  that  the  laminar  sublayer  has  a  thickness  proportional  to 
the  inverse  three  half  power  of  the  Reynolds  number  based  on  the  frictional 
velocity,  and  the  diameter  d. 

4.  The  Law  of  the  Wall  for  Injection  or  Suction 

The  shear  distribution  in  the  wall  region  with  injection  or  suction  is 
determined  for  a  flat  plate  with  zero  pressure  gradient  from  the  equation, 

.  du  1  dT 


wdy  p  dy 


Integrating  (34)  yields 


which  in  the  present  notation  may  be  written  as 


(ur)3  =vwu+  +  1 


where  v  =  v  /u 
w  w 

From  Eq.  (18) 
u 


r  du  '  r  du'  2  r,,  .  +  +A 

u  =  — r  =  — — i — r  -  — r-Ki+vu  )  2-i] 

i  dr4  i(vV+D%  v* 


and  the  integral  term  in  Eq.  (27)  becomes. 


4 


\ 


+ ,, 

_  P  du  P  du  l.r,.4-+1  2  .  Ti  .  W  T 

a "  J  — c  J~+T  “  "+6n[14vwu  3  =  r3 


(37) 


o  u_ 


w 


Substituting  Eq.  (36)  into  Eq.  (24)  yields  the  law  of  the  wall  for  injection  or 
suction  in  the  form, 


v  *  U 

y+  =  -^tnCl-H-j-]  +  |H[exp(klU)J  -  [l+(l^ U)+|(k1  U)2-h|(l^ U)3]|  (38) 

v 

w 

where  U  is  given  by  Eq.  (36).  Now,  either  from  the  asymptotic  form  of  Eq.  (38) 

or  directly  from  Eq.  (20),  the  region  outside  the  laminar  sublayer  and  the  buffer 

layer  is  described  by  the  equation 

U  =  —  ‘tnk8y+ 
ki 

Written  explicitly  in  terms  of  u+  as  determined  by  Eq.  (36),  EqJ  (20)  becomes 


•2-[Q4vV)“-1]  =  £  tnkay+  (35) 

v+  w  h 

w 

which  is  identically  the  modified  law  of  the  wall  as  suggested  by  Stevenson'*'0 
based  on  his  experimental  data. 

In  Fig.  (4)  the  law  of  the  wall  for  injection  or  suction  is  shown  for  four 

different  cases  (a)  v+  =  0.45  (b)  v+  =  0.064  (c)  v+  =  -0.051  and  (d)  v*"  ■  -0.066. 

WWW  w 

In  Figs.  (4a, b)  the  data  of  Stevenson  is  included  while  in  Figs.  (4c, d)  the  data 
from  Black  and  Sarnecki^,  as  measured  by  Black/  is  presented.  While  in  either 
the  large  injection  or  the  large  suction  cases  agreement  is  satisfactory,  in  the 
law  cases  where  experimental  scatter  seems  to  be  more  pronounced,  the 

agreement  is  at  most  fair. 

It  is  important  to  note  that  as  the  injection  rate  increases  the  buffer 
region  swells  on  the  lower  side  of  the  distribution  while  the  wake  component 
appears  at  lower  values  of  the  parameter  U, thus  leaving  an  extremely  narrow 
region  for  the  logarithmic  region.  It  seems  however  that  as  long  as  such  a 
region  exists  regardless  of  how  narrow  it  i3, the  buffer  layer  is  satisfactorily 
represented  by  the  present  theory. 
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SUMMARY 


A  generalized  law  of  the  wall  and  an  eddy  viscosity  model  for  wall 
boundary  layers  has  been  derived  for  regions  of  the  flow  field  which 
depend  only  on  a  single  space  variable.  The  derivation  is  based  on  Reichardt's 
result,  viz.,  mat  the  turbulent  shear  decays  as  the  cubic  power  of  the  distance 
from  the  wall  and  on  the  validity  of  the  logarithmic  law  of  the  wall  for  a 
constant  shear  flow  fx  Id. 

Considering  the  velocity  as  the  independent  variable  and  expressing 
€  as  a  function  thereof  yields,  in  the  constant  shear  case,  a  closed  form 
solution  for  the  velocity  distribution.  The  resulting  distribution  is 
shown  to  be  in  good  agreement  with  experimental  data.  Included  in  the 
comparison  are  some  recent  experiments  where  special  care  has  been  taken 
in  order  to  obtain  accurate  measurements  in  the  region  very  close  to  the 
wall. 

By  introducing  a  transformation,  an  extension  of  the  above  formulation 

to  variable  shear  has  been  obtained.  The  generalized  form  has  been  then 

applied  to  the  investigation  t  two  problems  (1)  the  effect  of  Reynolds 
* 

number  (u  d/v)  on  the  velocity  di  trillion  in  a  oipe  and,  (2)  the  deri  “ion 
of  the  law  of  the  wall  for  boundary  layers  with  injection  or  suction. 

In  the  first  problem,  it  has  been  demonstrated  that,  although  to 
a  very  small  extent  there  is  nonetheless  a  steepening  effect  the 
slope  of  the  logarithmic  portion  of  the  law  of  the  wall  with  an  increase 
in  the  Reynolds  number.  This  result  is  in  agreement  with  Hinze's  analysis 
of  Nikurudse's  data.  Furthermore,  an  expression  for  the  thickness  of  the 
laminar  sublayer  has  been  obtained  which  shows  that  (y  /d)  is  inversely 
proportional  to  the  3/2  power  of  the  Reynolds  number  (u  d/V). 

In  the  second  problem  a  law  of  the  wall  for  injection  or  suction 

+ 

has  been  derived  in  terms  of  the  parameter  v  .  Good  agreement  with 
available  experimental  data  has  been  indicated  and  a  tendency  for  a 
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decrease  in  Che  extent  of  the  logarithmic  portion  of  the  law  of  the  wall 
with  increase  of  injection  has  been  observed. 

From  the  experimental  data  considered, it  seems  that  the  model 
proposed  above  correctly  predicts  the  velocity  distribution  in  the  wall 
region  for  a  variable  shear  flow.  Obviously, a  more  complete  description 
requires  the  matching  of  these  results  with  an  eddy  viscosity  for  the  outer 
region. 
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PROBE  TIP  LOCATED  O.l mm. 


FIG.  ?  A  COMPARISON  BETWEEN  THE  PRESENT  ANALYSIS  AND  THE  EXPERIMENTS  OF  LINDGREN 


FIG.  3  VARIATION  OF  THE  EDDY  VISCOSITY  IN  THE  WALL  REGION  CLOSE  TO  THE  WALL  (ABBRECHT  AND  CHURCHILL 


v„*=+0.45(REF.  10) 


FIG.  4a  LAW  OF  THE  WALL  FOR  INJECTION  OR  SUCTION 


FIG.  4b  LAW  OF  THE  WALL  FOR  INJECTION  OR  SUCTION  (v  =  +0.064) 


FIG.  4d  LAW  OF  THE  WALL  FOR  INJECTION  OR  SUCTION 
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